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Abstract. We present, for the Ising model on the Cayley tree, some explicit formulae 
of the free energies (and entropies) according to boundary conditions (b.c). They in- 
clude translation-invariant, periodic, Dobrushin-like b.c, as well as those corresponding 
to (recently discovered) weakly periodic Gibbs states. The later are defined through a 
partition of the tree that induces a 4-edge-coloring. We compute the density of each 
color. 
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On non-amenable graphs, not only Gibbs measures but also the free energy (and the 
entropy) depend on the boundary conditions. 

The purpose of this paper is to study this dependence for one of the simplest such 
graph, the Cayley tree (Bethe lattice). Our analysis is restricted to the Ising model. 

Let r fc = (V, L) be the uniform Cayley tree, where each vertex has k + 1 neighbors 
with V being the set of vertices and L the set of edges. 

On this tree, there is a natural distance to be denoted d(x,y), being the number of 
nearest neighbor pairs of the minimal path between the vertices x and y (by path one 
means a collection of nearest neighbor pairs, two consecutive pairs sharing at least a 
given vertex). 

The Ising model is defined by the formal Hamiltonian 



where the sum runs over nearest neighbor vertices (x, y) and the spins a(x) take values 
in the set $ = {+1, -1}. 

For a fixed x° G V, the root, let 

W n = {x G V : d(x, x°) = n}, V n = {x G V : d(x, x°) < n} 

be respectively the sphere and the ball of radius n with center at x° , and for x G W n let 

S(x) = {yeW n+1 :d(y,x) = l}, 

be the set of direct successors of x. 



1. Introduction and definitions 




(1.1) 



(x,y)cV 
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The (finite-dimensional) Gibbs distributions at inverse temperature /3 = 1/T are de- 
fined by 

Mn(o-n) = Z' 1 exp |/3J ^ a(x)a(y) + ^2 h x a(x)}, (1.2) 
with partition functions given by 

Z n = Z n (p,h)= exp{/3J <?(x)cr(y)+ Y h * a ^}- ( L3 ) 

CTri G# v ™ (i,j)cv„ iew„ 

Here 

/i = {h x € i?, x G V} 

is a collection of real numbers that stands for (generalized) boundary condition. 
The probability distributions ( |1.2[ ) are said compatible if for all o n ^\ 

^ Mn(o"n-l,W n ) = /i n _i(cr„_i). (1.4) 

It is well known that this compatibility condition is satisfied if and only if for any x G V 
the following equation holds 

K= f( h v> e )> (1-5) 

where 

= tanh(/3J), /(/i, 0) = arctanh(0tanh/t). (1.6) 

Namely, for any boundary condition satisfying the functional equation (1.5) there 
exists a unique Gibbs measure, the correspondence being one-to-one. 

We will be interested in the dependence with respect to boundary conditions of the 
free energy defined as the limit 

F(P,h) = - lim * ]nZ n 09,/0, (1-7) 

TWOO fj\ V n \ 

where | • | denotes hereafter the cardinality of a set. 

A boundary condition satisfying (1.5) will be in the sequel called compatible. 

The paper is organized as follows. Section 2 provides the first part of results: a gen- 
eral formula applied then to various known boundary conditions (translation- invariant, 
Bleher-Ganikhodjaev, Zachary, ART), and those about entropy. Periodic and weakly 
periodic cases are the subject of Section 3. A first appendix concerns the density of 
colors mentioned in the abstract. A second appendix provides a sufficient condition for 
the existence of the free energy in case of compatible boundary conditions. 

2. General formula and first results 
For compatible boundary conditions, the free energy is given by the formula 

F(f3,h) = - lim -L V a(s), (2.1) 
n— >oo y„ * — » 

1 1 x£V n 
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where 

a(x) = — In [4 cosh (/i x - (3 J) cosh(/i x + pV)]. 
2p 

To see it, first notice that 

Z n (f3,h) = A n _ 1 Z n _ 1 (p,h), (2.2) 

where A n = F] b(x) with b(x) satisfying 

xew n 

] J exp(/3 Jew + uh y ) = b(x) exp(eh x ), e = ±1. (2.3) 

y&S(x) u=±l 

This formula used with both values of e implies 

6(a;)= ]"[ (4cosh(/ij / -/3J)cosh(/i^ + /3J)) 1/2 = exp(/3 ^ a(y)' 
ye5(x) yeS{x) 



It is then enough to insert this formula into the recursive equation (2.2 ) to get by iteration 

Z n ((3,h)= H b( X ) 



X&V n 



which gives (2.1 ). 
Notice that 

F((3,h) = F((3,-h), (2.4) 

where —h = {—h x ,x £ V}. 

2.1. Translation-invariant boundary conditions. They correspond to constant func- 
tions, h x = h, in which case the condition (|1.5|) reads 



h = kf(h,6). (2.5) 

The equation (2.5) has a unique solution h = 0, if 9 < 9 C = t an d three distinct solutions 
h = 0, ±h* (/i* > 0), when 6* > 9 C . 

Let us denote by no, \x± the corresponding Gibbs measures and recall the following 
known results for the ferromagnetic Ising model {9 > 0): 

(1) If 9 < 9 C , no is unique and extreme. 

(2) If 9 > C , and /^ + , are extreme. 

(3) /xo is extreme if and only if 9 < A=. 

(see e.g. |], §, [|) 

According to formula (2.1), the free energies of translation-invariant (TI) b.c. are 
given by: 

Fti(/3,0) = -iln(2cosh(/3J)). (2.6) 

FtiOSA) = F T1 (/3,-h*) = -— ln[4cosh(pV-/j*) coshOS J + /**)]. (2.7) 

zp 

Some particular plots are shown in Fig. [TJ 
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Figure 1. The free energies i*Tl(/3, 0) (solid line) and F T i(/3, /i*) (dotted 
lines) for >9 C with J = 1 and k = 4, 5, 6. 



In order to draw the free energy Fti{P,K) as a function of /3, we notice that the 
equation (2.5) gives 

a (i+|)2h 



1 



In 



1 



2 J e 2/i 



2h 

e k 



21 



and use the parametric representation defined by the following mapping: 

h -> (/3(/i), F(/i)), for /» > 0. 



The function F(h) is defined by inserting (2.8) into (2.6) and (2.7). 



2.2. Bleher-Ganikhodjaev construction. Here one consider the half tree. Namely 
the root x° has k nearest neighbours. Consider an infinite path tt = {x° = xo < x\ < . . . } 
(the notation x < y meaning that pathes from the root to y go through a;). Associate to 
this path a collection h w of numbers given by the condition 

f — h*, if x -< x n , x G W n , 

\ h*, if x n -< x, x G W n , 



hi 



(2.9) 
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n = 1, 2, . . . where x -< x n (resp. x n -< x) means that x is on the left (resp. right) from 
the path it. 

For any infinite path it, the collection of numbers h n satisfying relations (1.5) exists 
and is unique (see j2|). 

A real number t = t(n), < t < 1 can be assigned to the infinite path and the set 
is uniquely defined. The set of numbers h n ^' being distinct for different t £ [0, 1], 
it is also the case for the corresponding Gibbs measures. One thus obtains uncountable 
many Gibbs measures and they are extreme. 

Since the solution /i 7r can differ from or —h* only on the path tt, it is thus obvious 
that the values /ij, x S it do not contribute to the free energy. As a consequence, we get 
from the evenness (2.7), that the free energies of Bleher-Ganikhodjaev and translation- 
invariant boundary conditions coincide: 

F BG ((3,h*) = F Tl ((3,K). (2.10) 

2.3. Zachary construction. This con struc tion provides an (uncountable) set of dis- 
tinct functions h^> satisfying condition (1.5) and parameterized by t G (— /i*,/i*). It is 
assumed here that 9 > 9 C . 

Take t / and define the sequence (t n )n>o recursively by to = t, 

t n = kf(t n+1 ,9), n>0. (2.11) 

Since the function / is increasing and maps the interval (— /i*,/t*) into itself, the definition 
of t n make sense. Moreover one can see that lim n _ s>00 t n = for each to = t- 

Consider the function hx = t n for all x 6 W n . This function satisfies condition 



(1.5) for any t and by construction, distinct t assign distinct h^K The associated Gibbs 



measures are known to be extreme |11| . 

The corresponding free energies can be written as 

1 1 n 

*ZachCM (t) ) = -KB lim uT\ Yl \ W m\ °(tm] 
2p rwoo \ V n \ ' 
m=0 

where a(t) = ln[4cosh(/3J - t) cosh(/3J + t)}. 

By Stolz-Cesaro theorem (see e.g. [7]) applied to the sequences 



a n=^2 \ W ™\a(tm), b n = \V n \ (2.12) 



m=0 



one has 



n+1 



\W m \ a(t m ) - \ W m \ a(t m ) 

lim — - y \W m \ a(t m ) = lim 



n— >oo 



m=0 m=0 

\V n \ ^— ' ' "" V '"" «->oo \V n +l\ - \V n \ 

in=0 



hm l^+il "(Wi) = lim ~ a{tn+l) = ~ a{0) _ (2 . 13) 
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As a consequence 

F Zach (l3,hW) =F TI (/3,0). (2.14) 



2.4. ART construction. Let h be a boundary condition satisfying (1.5) on r feo . For 
k > ko + 1 define the following boundary condition on T k : 

{ h x , if x € V k ° 
h * = { n * ( 2 - 15 ) 

[ o, if x g 

where y fc denote the set of vertices of T k . Namely, to each vertices of V k ° one adds 
k — ko successors with vanishing value of the boundary condition. It is obvious the b.c. h 
satisfy the compatibility condition (1.5). In this way one constructs a new set of Gibbs 
measures that are extreme in the range 1/ko < 9 < 1/Vk [lj. 
For the corresponding free energy, we have 

F ART (P,h) = ~lim ^((|Kf|-|K fc0 |)ln[2cosh(/3J)]+/3 £ a(x)) . (2.16) 
Since 



x&V„ 



Ih^hllnq fo + l)*ff-2 

n->oo \V k \ k -l n^oo (k + l)k n - 2 

by taking into account < a(x) < C&, we get 

0< Yl a{x)<\V*°\C p . 



As a consequence, 

lim 

\v k \ 



lim Trm Yj a ( x ) = °> 



so that 

F A RT(P,h) = -^ln[2cosh(/3J)] =F T i(/3,0). (2.17) 

2.5. Entropy. To compute the entropy S(f3, h) = — dF j^' h ^ , we first notice that equation 
d2i8b gives 



1 cosh (2h) - cosh (fO 

h'{p) = —— = jk — y^. 2.i8 

y ' P'(h) sinh (2h - fcsinh f ) V ; 



As a result of easy computations, we get the formula 

S(fi, h) = -ln [2 cosh(2/i) + 2 cosh(2/3 J)] 



sinh(2h)-fcsinh(^) 



cosh(2/i) + cosh(2/3 J) 
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and 

S(P, 0) = ln(2 cosh(/3J)) - (3 J tanh(^ J) . (2.20) 



Let us mention that in case k = 2,we get by solving equation (2.5) the following value 
of /i* as a function of the inverse temperature: 



± h* = — In 



2- 1 (V^ - 2e 2/3J - 1 ± (e 2/3J - l)y/(e 2 P J + l)(e 2 ^ - 3) 



(2.21) 



The free energies and entropies then read 



F(J3,h*) = -i(ln(e 2 ^-l) + ^ln(e- 2 ^ + l)), (2.22) 

S(J3,K) = ln(2sinh(/3J)) + \ ln(2 cosh(/3 J)) - PJ *^^}* 1 ■ (2.23) 

2 2smh(2pJ) 

3. Periodic and weakly periodic Gibbs measures 

3.1. A group representation of the Cay ley tree. Let Gk be a free product of k + 1 
cyclic groups of the second order with generators ai, 02, . . . , Ofc+i, respectively. 

It is known that there exists a one-to-one correspondence between the set of vertices 
V of the Cay ley tree T k and the group Gk- 

To give this correspondence we fix an arbitrary element xq E V and let it correspond 
to the unit element e of the group Gk- Using a%,..., a^+i we numerate nearest-neighbors 
of element e, moving by positive direction (see Fig. [2]). Now we shall give numeration 
of the nearest-neighbors of each aj, i = 1, . . . , k + 1 by ajOj, j = 1, . . . , k-\- 1. Since all aj 
have the common neighbor e we give to it ajOj = a 2 = e. Other neighbors are numerated 
starting from by the positive direction. We numerate the set of all nearest-neighbors 
of each aiaj by words aiCija q , q = 1, . . . , k + 1, starting from a^a^aj = a, by the positive 
direction. Iterating this argument one gets a one-to-one correspondence between the set 
of vertices V of the Cayley tree T k and the group Gk- 

In the group Gk, let us consider the left (right) shift transformations defined as follows. 
For go £ Gk, let us set 

T g (h)=gh, (T g (h) = hg), for all h E G k . (3.1) 

The set of all left (right) shifts in Gk is isomorphic to the group Gk- 

3.2. Periodic boundary conditions. In this subsection, we consider periodic solutions 



of (1.5) and use the above group structure of the Cayley tree. 

Definition 1. Let G be a normal subgroup of the group Gk- The set h = {h x : x E Gk} 
is said to be G-periodic if h yx = h x for any x E Gk and y E G. 

Let 

,(2) 



G\ = {x E Gk ■ the length of word x is even}. 

k 



(2) 

Note that G\ is the set of even vertices (i.e. with even distance to the root). Consider 



the boundary conditions h and h T : 




, K, if x g Gf ] 
~K = { , , (3-2) 

if xeG k \Gf\ 



and denote by /J>^) the correspon ding Gibbs measures. 



The G- periodic solutions of equation ( 1.5 ) are either translation-invariant (G^-periodic) 
or G*jp -periodic (see J5J), they are solutions to 

u = kf(v,9), v = kf{u,9). (3.3) 
In the ferromagnetic case only translation invariant b.c. can be found. In the antifer- 



romagnetic case (9 < 0) the system (3.3) has a unique solution h = if 9 > —1/k, and 
three distinct solutions h = 0, and if 9 < —1/k. 

Let us also recall that for the antiferromagnetic Ising model: 

(1) If 9 > — 1/ k, is unique and extreme. 

(2) If 9 < — 1/k, and are extreme, 
see [6]. 

For periodic measures, we have 

FperCS./i^) = -~ ln[4cosh(/3J-/u)cosh(/3J + hj\ = F Tl (f3,K). 

3.3. Weakly periodic Gibbs measures. Let Gk/Gk = {Hi, H r } be a factor group, 
where Gk is a normal subgroup of index r > 1. 
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Definition 2. A set h = {h x ,x G G^} is called - weakly periodic, if h x = hu, for 

any x G H \ , xy G Hj , where x^ denotes the ancestor of x. 

Weakly periodic b.c. h coincide with periodic ones if h x is independent of x\_. 
Here, we will restrict ourself to the cases of index two and recall that any such subgroup 
has the form 



H A = < x G G k : ^ 
{ i&A 



even 



where = {1, 2, . . . , k + 1}, and u x {ai) is the number of ai in a word x G G k . 

We consider A ^ N^: when A = N k weak periodicity coincides with standard periodicity. 
Let Gk/HA = {Hq,Hi} be the factor group, where Hq = Ha,H\ = Gk \ Ha- Then, 



h x 



(3.4) 



(3.5) 



in view of (1.5), the .f/^-weakly periodic b.c. has the form 

/ii, x G H , x± G H , 
h 2 , x G H , xi G Hx, 
h 3 , x G Hi, xi G H , 
, hi, x £ Hx, xi G Hx, 
where the hi satisfy the following equations: 

hx = \A\f(h 3 , 9) + (k -\A\)f(hx,9), 
h 2 = (\A\ - l)f(h 3 , 0) + (k + l- \A\)f(hx,e), 
h 3 = (\A\ - l)f(h 2 , 9) + (k + 1 - \A\)f(h A , 9), 
h 4 = \A\f(h 2 ,9) + (k-\A\)f(hi,9). 
For sake of simplicity, consider k = 4 and \A\ = 1. In this case 

Hq = {x G Gk : LU x (ax) is even}, 

Hx = {x G Gfc : u) x (ax) is odd}. 

These sets are shown in Fig. [3| 

Let us recall the following results of [9J: 

There exists a critical value a CT (~ 0, 1569) of a = e~ 2 ^ J such that: 

(1) If a > a CT , there exists a unique weakly periodic state /io, 

(2) If a = a cr , there are three distinct weakly periodic states fio, fi^ , [if . 

(3) If < a < a cr , there are five distinct weakly periodic states [Xq, ^ > IH ■ 
These measures correspond to solutions of system (3.5) on the invariant set hx = —h±, 
h 2 = —h 3 , that is to solutions of: 



3f(hx,9)- f(h 2 ,9), 
4f(hx,9). 



(3.6) 



More results about weakly periodic Gibbs measures can be found in (9~1 , [To] . 
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Figure 3. The sets Hq (black vertices) and Hi (gray vertices). 

Denote 

Ai = \{{x,y) G L n : x G H , y = x± G H }\ , 
&n = \{(x,y) G L n : x G H , y = x± G H ± }\ , 
C n = \{{x,y) G L n : x G Hi, y = x± G # }| , 
£>n = G L n : x G Hi, y = x± G #i}| , 

where L n is the set of edges in V n . 
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For weakly periodic b.c. (3.6) we have 

Fyrp(p,h) = -— lim -^{(^ n +^n)ln[4cosh(/3J-/ii)cosh(/3J + /ii)] + 

2p n->oo | V n \ 

(B n + C n ) ln[4 cosh(/3 J - h 2 ) cosh(/3 J + h 2 )]} ■ 
By Proposition [T] of the appendix we obtain 

F WP (p,h) = (^ln[4cosh(J/3-/ti)cosh(J y 9 + /t 1 )] 

+ ~ ln[4 cosh( J/3 - h 2 ) cosh( J/3 + /i 2 )] j (3.8) 

In the case under consideration the proposition is straightforward. Indeed one imme- 
diately see, in view of definitions and Fig. [3| that (A n + V n )/(\V n \ - 1) = 4/5 and 
{B n + C n )/{\V ri \ — 1) = 1/5 when n = 1; the induction is trivial. 
The equation 

hi = 3f(hi,0)-f(4f(h 1 ,O),9) (3.9) 
that solves the system ( |3.6| ) (with /i 2 = 4f(h\,9)) can then be reduced to: 

a 2 f - af - 2o?i + a + 1 = 0, (3.10) 

which has two solutions £i and £2 when < a < a cv (see [9]). 
The free energy then reads: 



( 



1 

10/3 \ 



4 In 



w „ « 2^ cosh 2/3 J) - 4 

2cosh(2/3J + / „J — - 

v P ; 2cosh(2/3J -£ 



+ ln[2 cosh(2/3 J) + £ 4 - if + 2] }. 



The solution £1 is given by 

1 / 2ft 2 1 /3( 6 + e 4/3) \ 

a = 3 1 e 2/3 + — — + 2 ^ J ' 

where 

U = (-9e 2/3 - 27e 4/3 + 2e 6/3 + 3^-96 - 39e 4 ^ + 54e 6 " + 69c 8 / 3 - 12e^) 1/3 . 
The corresponding free energy reads 

1 (_ e -2/3 + e 2pf (2 + e" 2 ^ + e 2 / 3 - 4 (e 4 ^) 2 + (e 4 ^f) 4 ) 
^WP(^) = -^In (e - 2 , + e2 ,_ e4 , f)4 ' 

where 

W = (2e- 2 ^ + ^ + 2 2 / 3 e" 4 %) with V = 2^e~^ (6 + e 4 ' 3 ) . 

The plot is given in Fig. [4] from which we observe the strict inequalities 

F T i(/3, K) < F W p(/3, fc) < F T i(/3, 0). (3.11) 
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in the range a < a CI 
o 




Figure 4. The free energies i*wp(/3) (dotted line) for a < a CT to- 
gether with the previous free energies Fti(/3, 0) (solid line), and -Fti(/3, h*) 
(dashed line). Here J = 1 and k = 4. 

The results for weakly periodic boundary conditions have to be compared with the 
inequalities given recently in [4j; the weakly periodic b.c. (3.6) corresponding to the 
so-called dimer covering in (4j. There, the inequalities are easier to catch with cluster 
expansion method in mind, the condition on the temperature is more restrictive, and 
the free energies cannot be expressed explicitly. 

Appendix: Density of edges in a ball 

In this appendix we consider a group representation of a Cayley tree and its partition 
with respect to an arbitrary subgroup of index two. This partition gives a 2-vertex- 
coloring on Cayley tree, which then gives 4-edge-coloring, say, colors i = 1,2,3,4. We 
fix a root of the Cayley tree and give explicit formulas for number A n ^i of edges with 
color i in a ball V n of radius n with the center at the root. Moreover, we compute the 

lim,i_ K30 (.A J i )i /|K|) for each i = i' 2 ' 3 ) 4 - 
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We will use the notation of Subsection 13.31 and let 

a n = \{{x, y) £ L n \ L n _i : x £ H , y = x i £ H }\ . 

Pn = \{{x,y) £ L n \ L n _i : x £ H , y = x± £ i?i}| . 
In = \{{x,y) £ L n \ L n _i : x £ H 1 , y = xi £ H }\ . 
5 n = \{{x,y) £ L n \ L n „i : x £ Hi, y = x± £ Hi}\ , 

for A = {l,2,3,...,j} and 1 < j < k+1. 

Let M be the set of all unit balls with vertices in V and let S\(x) denotes the set of 
all nearest neighbors of x. For b £ M the center of b is denoted by c b . 

Lemma 1. If Cb £ Hq, then 

\{x £ Si(c b ) : x £ Hi}\ = j, \{x £ Si(c b ) : x £ H }\ = k - j + 1. 

// Cb £ Hi, then 

\{x £ Si(c b ) : x £ Hi}\ = k - j + 1, \{x £ Si(c b ) : x £ H }\ = j. 

Proof. We have Si(c b ) = {c b a p : p = 1, 2, . . . , k + 1}. 

If c b £ Ho, (the case c b £ Hi is similar) then c b a p £ Hi for p = 1,2, and c b a p £ Ho 
for p = j + 1, j + 2, k + 1, i.e. we have 



\{x £ Si(c b ) : x £ Hi}\ = j, \{x £ Si(c b ) : x £ H }\ = k - j + 1. 



□ 



(A.l) 



Consider b = V\ £ M with the center x° = e £ Hq, then in Wi we have j vertices 
which belong to Hi, and k — j + 1 vertices which belong in Hq, consequently, 

ai = k-j + l, (3i = 0, 7i =j, 5i = 0. 

Lemma 2. For any n £ N the following recurrence system hold 

a n+ i = (k- j)a n + (k- j + l)/3 n 

Pn+l = (j ~ l)7n + j^n 

7n+l = (j - l)/3 n + ja n 
k <5„+i = (k - j)5 n + (k- j + 1)7„, 
with initial values ai = k — j + 1, Pi = 0, 71 = j, 5i = 0. 

Proof. By Lemma [TJ an edge (x,y) £ L n \ L n -i with x £ Hq, y = x^ £ Hq has (A; — j) 
neighbor edges (z, x) £ L n+ i \ L n with z £ Hq, x = z± £ Ho. An edge (z, t) £ L n \ L n _i 
with z £ Hq, t = zi £ Hi has (k — j + 1) neighbor edges (u, z) £ L n+ i \ L n with 
u £ Hq, z = £ Hq. Moreover, it is easy to see that only a n and /3 n have contribution 
to a n +i. Hence we have a n +i = (k — j)a n + {k — j + l)p n . Other equations of the system 
(A.l) can be obtained by a similar way. □ 



In 



, m = l,2,... 



, m=l,2, ... 
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Remark 1. For j = k + 1 by Lemmas^ and[S] we get a n = 5 n = 0, for any n > 1 and 

{0, if n = 2m — 1 
(/c + l)*: 2 " 1 - 1 , t/ n = 2m 

0, if n = 2m 

(k + l)^" 1 -^, if n = 2m-l 
So in the sequel of this section we consider j as 1 < j < k. 

Lemma 3. For a n we have 

a n+2 = j(k -j + l)\W n \ + (k — 2j)a n+1 - (2j - l)a n - ka n -i, n>2, (A.2) 

with initial values 

ai = k-j + l, a2 = (k-j)(k-j + l), a 3 = ((k-l) 2 +j(j-l))(k-j + l) (A.3) 

Proof. The initial values follow from Lemma [2} 
By definitions of a n , (3 n , j n , 5 n we have 

a n + Bn+ln + S n = \W n \=k n - 1 (k + l), U > 1. (A.4) 



From (A.l) we get 

Pn = k^j+i{a n+1 -{k- j)a n ), 

In = k-jA-1 {o-n ~ (k-j)an-l) + ja n -l, 

5 n = \(p n+l -{j-l) ln ). 
Substituting these values in ( |A,4 ) and then simplifying we get (A.2). 



(A.5) 



□ 



To find solution of (A.2) we denote 



n-2 



(A.6) 



From (A.2) we get 

k n q n+2 = k n -\k +l)(k-j + l)j + (k — 2j)k n ~ 1 q n+1 - (2j - l)k n - 2 q n - k n - 2 q n ^, 



dividing by k n we obtain 



q n +2 = 
with initial values 



(k + l)(k-j + l)j , k-2j 



k 



+ 



k 



ln+1 



2j - 1 J_ 

k* Qn 



ln-1, 



Ql = + Q2 = (k - j)(k - j + 1), q 3 



((fc-l) 2 +j(j-l)) (k-j + 1) 



(k+i)(k-j+i)j 



In order to find solution to (A.7) first we rid j- 

, Kk-j + i) 

q n =Pn-\ 7, 



A: 

by denoting 



(A.7) 



. (A.8) 



(A.9) 



ON FREE ENERGIES OF ISING MODEL 



15 



Substituting (A. 9) to (A.7) we get 



Pn+2 



k-2j 



j 2j-l 1 



(A.10) 



with 

1 (fc-2j)(fe- j + 1) (fc 2 -4fc + 2 + 2j 2 -2j)(fc-j + l) 
Pi = ~3 + V,P2 = ^ , P3 = ^ • 

(A.ll) 



The characteristic equation for (A. 10) has the following form (setting p n = A™): 
which has solutions 



, 1 , k - 2j + 1 ± y/(k - W - 2(fc + 2jJTl 

Ai = --, A 2i3 = ^ . (A.12) 



Then the general solution to (A. 10) is 



Pn = AiA" + A 2 A£ + A 3 Xl 



(A.13) 



where the coefficients A\, A 2 , A3 are determined by the initial conditions (A.ll). 



Using ( A.9 ) and ( A. 6 ) we get 



(-1)™ , M 



(A.14) 



Then using (A. 5) and (A.14) one can find /3 n , 7„ and <5 n . 
We have 



m=l 

A 2 A 2 £; 



«w - 1) + ,^ 3 y ^ (( w - 1) • 



k 2 (X 2 k- 1) vv z y ' 1 k 2 (X 3 k- 1) 
Proposition 1. For any j = 1, . . . ,k and any fixed q = 0, 1, 2, . . . we /iawe 
li m ^!!zi = i im 5 Jlzl = (fc- l)(fc- J + 1) 

lim = Iim 7 ^ - ( * " 1)j 



(A.15) 



n-+oo \V n \ n-+oo \V n \ 2(k + l)k1+ 1 ' 

k-j+l 

\V n \ 2(k + l)kfl' 



lim 

n— >oo 



lim 



A n — q 


= lim 




\Vn\ 


n— >oo 


\Vn\ 


$n—q 


= lim 


Cn—q 


\Vn\ 


n— yoo 


\Vn\ 



3 
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Proof. It is easy to check that IA2I < 1 and | A3 1 < 1, i.e. 



k-2j + l± ^{k - 2j) 2 - 2(fc + 2j) + 1 



< 2k, for any 1 < j < k. 



Using these inequalities, formula (A. 14) and (A. 16) we get 

(k-l)(k-j + l) 



a 



lim 

n— >oo v r 



n—q 



2(k + l)ki+ 1 



Now using this formula together with (A.5) we obtain 

fin—q _ 



lim 



1 



n— >oo \V n \ 



k-j + l 



lim 

n— >oa 



IKI 



Q: 



n—q 



j) lim 

n— >oa V n 



J(fc-l) 
2(fc + l)^ 1 ' 



The formulae involving 7 n and S n are obtained in a similar way. 

(k + 1) • k n - 2 



By dAA5l) and 



IKI 



fc-1 



(A.16) 



we get 



1 . An— i 

hm —— 

71— !>00 \V n \ 



k-j + l 

2(k + l)kfl' 



From (A.5) we get 



k-j+l (Ai - «i + a n+1 - (k - j)An) 
jrrj^x (An ~ (k - j)An-l) + j'Ai-1, 

: i (#„ - /3i + /3 n+ i - (j - 1)C) , 



that allows to prove the remaining formulae. 



(A.17) 



□ 



Remark 2. By Proposition^ it is clear that the values of A\, A 2 , A3 do not give any 
contribution to the equalities of the proposition. This is why we did not compute A\, 
A2, A3. But one can obtain the numbers by the initial conditions (A. 11) for p n . For 

: = 4 

Pn = Ai(- 



example, in the case k = 4, j = 1 from (A. 13) and (A. 11) we have 



pi 



The initial conditions give 



-f + A 2 *=0 + A 3 

A k + A2 l=gji + A3 l+^7i 



16 



64 1 " i2 128 3 128 

Consequently, (for k = 4, j = 1) we have 




_l_ A -9-5V7i , a -9+5V7: 

+ 198 no 



20y/7i 
7 

20y/7i 
7 



a, 



2 • 4™" 1 + 



1 



1 




(A.18) 
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Using (A. 5) and |A. ltfy one can find (3 n , j n and 5 n . Moreover, we have 
, A _ 2(4" - 1) 2V7i (3 - V7i) n - (3 + y/7i)' 



3 7 2™ 

m=l 



Note that a n and A n are natural numbers for any n > 1 . 
Remark 3. In the case j = k + 1 by Remark [7] we get 



lim —— = lim ■— = lim — — - = nm 



\V n \ n— >oo \ V n \ n— s>oo |V^| ra— >oo 

,. (fam-X ,. 72m n 

hm TT7 1 = lim TT7 — T = °- 

m->oo |V2m-l| m ^°°|l / 2m| 

1- ^2m _ 72m- 1 _ k-1 

m->oo \ V2m\ m-+oc |V2m-l| & 

#2m C2m-X k 

lim — — r = lim 



m->oo |V2m| m-s-oo |V2m-l| + 1 

,. 02m-l t ^2m 1 

lim r = lim 



m-»oo V2m— 1 m-»oc I Vam k + 1 

Appendix: Existence of the free energy 

As we have seen in the previous sections free energy exists for each known compatible 
boundary condition. We note that a(x) is bounded: /3 _1 ln2 < a(x) < Cr. Hence limit 
(2.1) is also bounded. But the problem of convergence of (|2.1|) is still open. Here we 



shall give some conditions on h, under which the corresponding free energy exists. 

Let 7T = {x° = xq < x\ < . . . } be an infinite path. A function h x on the path 
7r is called monotone non increasing (non decreasing) if h x . > h Xi+1 , (h Xi < h Xi+1 ), 
» = 0, 1, 2, 

Proposition 2. Let h = {h x , x £ V} be a compatible b.c. If on any infinite path starting 
at x £ W no , \h x \ is monotone non increasing (non decreasing), then the corresponding 
free energy F(j3, h) exists. 

Proof. Using Stolz-Cesaro theorem we get 



lim — r V" a(x) = lim A n , with A n = — — r V" a(a 

1 1 xeV n 1 1 X€Wn 

We shall show that A n is monotone for n > uq. We have 

A„_i - A n = y— [ k E - E a ( x ) 



— E *»(*)- E \=w-\ H E K*) -«(*))■ 



18 



D. GANDOLFO, M.M. RAKHMATULLAEV, U. A. ROZIKOV, J. RUIZ 



By monotonicity of \h x \ and by evenness of a(x) we notice that a(x) — a(y) does not 
change sign for all x, y with x < y. Thus A n is monotone and since it is a bounded 
sequence it has a limit. □ 
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